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^3 : ABSTRACT 

1 1 A phenomenological formalism is presented in which the apparent acceleration 

,— il of the universe is generated by large-scale structure formation, thus eliminating 

O ■ the magnitude and coincidence fine-tuning problems of the Cosmological Con- 

stant in the Concordance Model, as well as potential instability issues with dy- 
namical Dark Energy. The observed acceleration results from the combined effect 
^ I of innumerable local perturbations due to individually virializing systems, over- 

^ I lapping together in a smoothly-inhomogeneous adjustment of the FRW metric, in 

^ I a process governed by the causal flow of inhomogeneity information outward from 

each clumped system. After explaining why arguments from the literature claim- 
ed ■ ing to place restrictive limits upon backreaction are not applicable in a physically 
ly-^ I realistic cosmological analysis, we present a selection of simply-parameterized 
^ ■ models which are capable of fitting the luminosity distance data from Type la 
if^ ■ supernovae essentially as well as the best-fit flat ACDM model, without resort to 
G^ ■ Dark Energy, any modification to gravity, or a local void. Simultaneously, these 
2 I models can reproduce measured cosmological parameters such as the age of the 
^ I universe, the matter density required for spatial flatness, the present-day decel- 
^ ■ eration parameter, and the angular scale of the Cosmic Microwave Background 
^ ■ to within a reasonable proximity of their Concordance values. A potential ob- 
■ servational signature for distinguishing this cosmological formalism from ACDM 
may be a cosmic jerk parameter significantly in excess of unity. 

Subject headings: cosmological parameters — cosmology: theory — dark energy 
— large-scale structure of universe 



1. INTRODUCTION AND MOTIVATIONS 



^Evidence from t he past decade or so indicating that we liv e in a spatially fiat ( iLarson et al, 

201ll ). accelerating ( iPerlmutter et al.lll999t iRiess et al.lll998[ ) universe, yet one in which the 
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density of clustering matter is significantly less than the critical density (e.g., lTurnerll2002l ). 
has led cosmologists to a surprising conclusion: the dominant cosmic component appears to 
be "Dark Energy", a smoothly-distributed ma terial that must posse ss negative pressure in 



order to be able to drive the acceleration (e.g., iKolb fc Turnerlll990l ) 



The true nature of this Dark Energy remains very much in doubt, however. The simplest 
version, a Cosmological Constant (A), suffers from well-known aesthetic difficulties, including 
two different fine-tuning problems: one being the issue that p\ is ^120 orders of magn itude 
smaller than what would be expected from the I ^lanck scale fe.g..lKolb fc Turner! 1 19901) : and 



the other being a "Coincidence Problem" (e.g., lArkani-Hamed et al. 



200ol ) which questions 



why, given pa/pm oc a^, that observers today should just happen to live right around the 
time (within a factor of ~2 in scale factor) when A began dominating the cosmic evolution. 

One possible way of ameliorating such problems is to give the Dark Energy an evolving 
equation of state, w(z), where such co incidences are eliminated or moderated (like in tracker 
quintessence; e.g.. IZlatev et al.l ( 119991 )). But moving away from a Cosmological Constant by 
giving the Dark Energy an evolution in w{z) also opens the door to other possible dynamics 
- such as the ab ility for the Dark En ergy to be mobile and thus perhaps to become spatially 
inhomogeneous (ICaldwell et al.lll998[ ). A dynamical Dark Energy (DDE) of this type is liable 
to join in structure formation to an observationally unacceptable degree, thus ruining it as 
a candidate to serve as the "smoothly-distributed" missing cosmic ingredient. 

From the 1^* Law of Thermodynamics, Pressure = — 9 [Energy] /9[ Volume], so that a 
'negative pressure' (P < 0) substance like Dark Energy will have dE/dV > 0; consequently, 
it should be self-attractive on local scales (and relativistically so, since w = P / p —1), in 
stark contrast to how Dark E nergy is popularly described as a "repulsive" or "antigravity" 
material (e.g., iHeavensI 120101 ). This potential instability of DDE to spatial perturbations 
must be countered through the addition of some nonadiabatic support pressure (an adiabatic 
positive p ressure term would cause cosm ic deceleration)] an d though this can app arently be 



done (e.g.. lCaldwell et al.lll998l : lHull2005[ ). other results exist (IMaor k. Lahavll2005l ) indicating 



that a Dark Energy allowed to dynamically participate in virialization does indeed show signs 
of providing an additional attractive force for clustering, rather than opposing clustering or 
leaving it unaffected. In any case, an inhomogeneous DDE supported against extensive 
clumping by nonadiabatic pressure would not be a perfect fluid, and since a perfect fluid 
approximation is th e basis of the argument for acceleration (a > 0) from P < (e.g., 
Kolb fc Turner! 1 1990l pp. 48-50), the entire situation is best avoided if possible. 



Interestingly, the self-attractive nature of negative pressure presents us with the seed 
of an alternative idea: in its own way, normal gravitational attraction represents a form of 
'negative pressure' - in particular, the gravitational attraction responsible for the growth of 
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clustering in large-scale structure formation. Making a virtue of necessity, one may there- 
fore recruit cosmic structure formation as the driver (not just the trigger, as in tracker 
quintessence) of the cosmic acceleration. This would completely solve the Coincidence Prob- 
lem, since the remarkably contemporary onsets of both the acceleration and the existence 
of galaxies (and hence planets, life, and astronomical observers) would be nothing more 
coincidental than finding two neighboring apples that have fallen from the same tree. 

Such considerations are what initially led this autho r to search for a structure-formation- 
induced solution to the cosmic acceleration problem (IBochnerl 120071 ) , though a plausible 
acceleration mechanism was not yet apparent. During the pas t few years, a number of other 



researchers have similarly propos ed a variety of methods (e.g.. iSchwarzl 120021 : lRasanerul2004l : 
Kolb et al.ll2005l : IWiltshird 120071 . as a few examples) also attempting to obtain the observed 
cosmic acceleration from the breaking of FLRW homogeneity and the formation of structure, 
without requiring the intervention of non-standard gravitational or particle physics, or the 
non-Copernican approach of placing us near the center of a large cosmic void. This type 
of approach has become generally known as "backreaction" . The backreaction paradigm, in 
the consensus view, has been unable to substantially replace Dark Energy as the source of 
the apparent acceleration; a conclusion which we will argue is quite premature. 

In this paper, in fact, we will present a physically plausible (though phenomenological) 
mechanism that is not only able to reproduce the cosmic acceleration as indicated by ob- 
servations of standard candles such as Type la supernovae, but which can also reproduce 
several other key features of the conventional ACDM concordance. 

The paper will be organized as follows: in Section [21 we explain why several oft-quoted 
limitations of the backreaction mechanism are inapplicable in realistic situations, and then 
use such arguments as guideposts towards developing our formalism in Sectional in Section HI 
we describe our specific models and present our simulated Hubble curves; in Section El we 
evaluate the success of our models in formulating an alternative "Cosmic Concordance" ; and 
we conclude in Section |6l Note that the in-depth technical details of the re search iiitrodu ced 
here are discussed much more fully in a companion paper by this author (IBochnerl |20 111 ). 



THE COSMOLOGICAL IMPACT OF "BACKREACTION" 



At present, backreaction i s not widely c onsidered to be a viable method for achieving 



the observed acceleration (e.g.. ISchwarzll2010l ). There are a number of important reasons for 



this view, reasons which seem on the surface to be strong arguments; but we will explain 
here why despite being technically correct, several of those arguments do not properly apply 
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to a physically realistic cosmological model. 

One popular approach for estimating (and placing limits upon) the effects of backreac- 
tion is to use "Swiss-Cheese" solutions - regular FLRW cosmologies with spherical 'holes' cut 
out of them, with precisely designed boundary conditions so that the exterior universe (the 
'cheese') remains completely unaffected by the holes. Such models are advantageous because 
exact solutions may be used for the clumped material inside the holes, allowing for a precise 
analysis. But the weakness of such models is their significant lack of realism: not only do 
most of them lack any mechanism for representing structure virialization and stabilization, 
but they have in fact been specially designed to set the actual backreaction (as opposed to 
purely observational effects) to zero, by hermetically sealing off each condensed structure so 
that it cannot affect the exterior universe (even shielding the holes from one another). In the 
real universe, though, the 'spheres of influence' of individual mass overdensities continually 
expand over time, affecting more and more distant objects, until multiple 'independent' influ- 
ences eventually overlap and merge. This superposition of gravitational perturbations from 
many different inhomogeneities is a phenomenon that Swiss-Cheese models cannot represent 
at all. Thus we understand why Swiss-Cheese backreaction calculations fail to find effects 
strong enough to produ ce the observed acceleration - such as, for example, the findings of 



Biswas &: Notaril (j2008[ ) showing that the perturbative effects of a hole on observations are 
strongly suppressed, by a factor {L/RhY for passage through the hole, where L is the size of 
the hole and Rh is the Hubble radius. Recognizing the improperly artificial act of abruptly 
cutting off the gravitational influence of an overdensity at hole 'edge' r = L, as the Swiss- 
Cheese models do, we see therefore that the smallness of the suppressive term [L/RhY is 
itself completely artificial. A realistic backreaction model would not have any predetermined 
limits on L, but rather would allow it to grow without bound (within causal speed limits 
for the propagation of gravitational perturbation information), overlapping with more and 
more distant holes, until the factor {L/RhY is no longer much of a 'suppression', after all. 

Another a rgument claiming to li mit the efficacy of backreaction is the acc eleration "no- 



go" t heorem of lHirata fc SeljakI ( 120051 ) . based upon the Raychaudhuri equation ( iHawking fc Ellis 
19731 ). demonstrating that a 'true' (volumetric) acceleration of a spatial region is impossible 
in the absence of a material with sufficiently negative pressure {w < —1/3) to serve as Dark 
Energy. This conclusion only holds, however, as long as one assumes the cosmic matter to be 
irrotational - that is, possessing zero vorticity. But here we must ask: Why is the vorticity 
set to zero? This seems to be an especially unusual prescription in a universe where nearly 
everything is rotating or revolving about something. Essentially all virialized structures in 
the universe larger than individual solid objects are stabilized against collapse by some ver- 
sion of vorticity or velocity dispersion; and it is certainly a bad physical approximation to 
take the dominant physical force (in opposition to gravity) which regulates cosmic structure 
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formation, and completely set it to zero. Sin ce vorticity in signific ant amounts - on large 
scales, at least - is not a predict ion of inflation (IHirata &: Seljakll2005[ ) , it is typically regarded 



as a "small-scale player" (e.g.. iBuchertl |2008[ ) . relevant only for cosmic averages performed 



over domains that are on or below the scales of gala xy clusters; and thus it is usually set 



to zero for technical convenience (e.g.. lRasanerul2010l ). a seemingly harmless simpliflcation. 



But this ignores the fact that the term actually appearing in the Raychaudhuri equation is 
not the vorticity tensor u,,,„, but th e squa re of the vorticity, w^, a quantity that will never 



average away. As iBuchert fc Ehlerd (Il997l ) put it, cosmic averages of positive semi-definite 
quantities like (cu^) get "frozen" at the size of "typical subdomains"; and thus they cannot 
be made to go to zero by averaging over larger domains, even if the spatially-averaged value 
of the parent tensor, {uj^^), does itself become negligible on large scales. Thus dropping 
vorticity entirely from the cosmic evolution is a very poor approximation. 

We do not attempt here to explicitly model the enormously complex backreaction effects 
of vorticity during the dynamical virialization of a structure; rather, we exploit the very 
simple natures of both the beginning state (nearly perfectly smooth FRW matter), and the 
ending state (reasonably randomly distributed Newtonian mass concentrations) of structure 
formation. This will lead (in Section [31 below) to a phenomenological formalism depending 
only upon the general evolution rate of the cosmic matter from "smoothness" to "dumpiness" 
over time. 

Lastly, note that our formalism uses the backreaction effects of Newtonian-strength (i.e., 
gravitatio nally linear) perturbatio ns, despite the extensive work of Buchert and collabora- 
tors (e.g.. iBuchert fc EhlersI 119971 ) seemingly demonstrating that the entire Newtonian-level 
backreaction can be expressed mathematically as a total divergence, thus guaranteeing it 
to be either extremely small or ider itically zero (depending u pon the boundary conditions). 
But cornparin g Equations la-d of IBuchert fc EhlersI (119971 ) to Maxwell's Equations (e.g., 
Jacksonlll975l ) shows that their formalism contains a crucial simplification: a dropping of all 
"gravitomagnetic" fields from their calculations. Similar to how structure formation calcula- 
tions are based upon the Poisson equation, = AnGp (though suitably modified for the 
cosmic expansion), in which the time-derivative terms of the full wave equation have been 
dropped, t his is an explicitly non- ca usal formalism. (Other backreaction no-go arguments 



like that in llshibashi &: Waldl ( 120061 ) are also reliant upon dropping terms like \d^/dt\'^, as 



seen via their Equation 2b.) This simplification is typically justified by citing the smallness 
of f^/c^ given the (reasonable) assumption of nonrelativistic speeds for most matter fiows; 
but there is a crucial difference between "small" and "suppressed" . Many individually small 
contributions can be added together to produce a large overall result, if one has enough of 
them; and though Newtonian-level perturbations weaken as ~l/r, the number of them in a 
spherical shell increases as (easily overwhelming any factors of f^/c^), creating a total 
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perturbative effect that would formally be infinite when integrated out to r = oo, if not 
reined in by the finite causal horizon out to which an observer can 'see' clumped structure 
that has had sufficient ti me since the Big Bang to form; a situation reminiscent of Olbers' 
Paradox (IWeinberg||l972[ ). 



For a causality-respecting approach, one must instead (as in electrodynamics) use the 
full wave equation, for special-relativistically consistent perturbation potential function ^sr- 

V2<l>SR-^^^ = 47rGp. (1) 

(Note that factors relating to the cosmic expansion are still neglected here, for simplicity.) 
Now, the usual impulse is to immediately drop the extra term in Equation [H involving 
d'^^sn/dt^ ^ equivalent to dropping the gravitomagnetic terms, as is done in the Buchert 
formalism - because of its resulting prefactor of v^/c^; this factor would seem to make it 
very small given the assumption of nonrelativistic speeds for most matter flows, and thus 
(assumedly) ensuring it to be negligible compared to the spatial variations term in any 
backreaction calculation. But this thinking is based only upon considerations of individual 
Fourier perturbation modes, not on the overall causal behavior of information flow in the 
structure-forming unive rse. If we inst ead keep all terms, and solve Equation [T] as-is, then 
one gets (adapting from I JacksonI llQTsl . eq. 6.69): 

$SR(x,t) = -G TM^^^dV, (2) 
J I-X -X I 

where the bracketed numerator is always evaluated at the retarded time, = t — |x — x'|/c. 
It is this retarded-time condition which restores causality, allowing different regions of the 
universe to communicate with (and gravitationally perturb) one another; and which provides 
the escape route from the Buchert suppression of Newtonian-level backreaction, because such 
backreaction is not truly expressible as a total divergence. We will refer to this propagation 
of gravitational perturbation information between distant (though communicating) regions 
as "causal updating". 

Given the fact that the key metric perturbation function, $sr(^); is predominantly af- 
fected by inhomogeneity information coming in from distant locations, the retarded-time 
condition of an integrated formula like Equation [2] (suitably modified for cosmological calcu- 
lations) implicitly gives it the ability to exhibit relativistic behavior in what would otherwise 
appear to be nonrelativistic situations. Consider the farthest distance (at some given obser- 
vation time, t) out to which the look-back time still represents a fairly clumpy universe; this 
distance (as a function of t) represents the edge of an outgoing "wave of observed dumpi- 
ness" expanding outward from an observer at c, which ensures that the potential function 
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$SR for perturbations felt by that observer will be fully (special-)relativistic in character, 
even if none of the matter in the universe is actually moving at relativistic speed^ This 
unexpected result is due to the fact that cosmological systems (unlike anything else) are 
infinite in size, meaning that the inherently relativistic act of causal observation will end 
up encompassing spatial volumes which are incomparably vast, causing even small effects 
to sum together to produce a dominant overall influence. Perhaps surprisingly, therefore, 
the apparent acceleration of the universe which we observe is not caused by any individually 
powerful or gravitationally exceptional specific objects or regions, but rather is an effect that 
comes from the summed influences of the weak-gravity Newtonian tails of innumerable mass 
concentrations, imposing their combined effects upon us from extraordinarily large distances. 

In this approach, it is clear that we are advocating a diametrically- opposed vi ew from 



that of other authors who have sought to define a "finite infinity" (e.g., iCoxl 120071 ). repre- 
senting a very limited boundary from within which significant effects upon some specified 
local volume can have ever arrived. Such a limitation would restrict the history of 'impor- 
tant' interactions to be within an effective "matter horizon" (lEllis fc Stoegerl 120091 ) of only 
a few Mpc in size, delineated by the timelike world line s trave led by pressure- free r natter 
due to scalar perturbations. But interestingly, both ICoxl (120071 ) and lEUis fc Stoegerl (120091 ) 
acknowledge the Great Attractor/Shapley Concentration as a likely source of our bulk flow 
- i.e., as the completely dominant influence upon the motion of our entire neighborhood of 
galaxy clusters - despite the clearly contradictory fact (for their claims) that those structures 
which control our bulk motion are exceedingly far beyond any reasonable estimate of our 
so-cal led matter horizon. Thus one cannot take at face value the assertion in lEUis fc Stoeger 
( I2OO9I ) that, "the strengths of any other possible long distance influences... gravitational 



waves, or electric Weyl tensor components from sources outsid e that region - ar e insignifi- 



cant compared to local effects." Notably, it has been shown by iLudvigsenI ( 1l989l ) that even 
an arbitrarily small energy flux due to gravitational waves can result in a finite amount 
of "geodesic deviation" (actually, long-term positional displacements of observers due to 
permanent metric changes) in the very distant radiation zone. 

Nevertheless, there are very straightforward calculations using gravitational perturba- 
tion theory which are commonly believed to contradict any assertions of the significance 
of such nonlocally-acting causal backreaction effects. The problem, however, is the fact 
that perturbation theory is singularly ill-equipped for dealing with the most important per- 
turbative effects in an unbounded system like the (effectively) infinite universe. The central 



^ As an analogy, consider the contact point between the two blades of a very long pair of scissors. The rate 
at which this contact point moves outward from the central pivot does not represent the physical motion of 
any real object, and hence is not limited by the speed of the material in the blades as they come together. 
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difficulty with perturbation theory here relates to its basic program of identifying the 'impor- 
tant' physics by labeling the amplitude of each term as "large" or "sm all" , and then drop ping 



the small terms in order to focus upon the large ones. For example, iKolb et al.l (|2005[ ) ne- 
glect information- carrying tensor modes generated by nonlinear scalar perturbations (i.e., 
virializing masses) in the final expressions for their anal ysis (as w e ll as a ssum ing irrot a tiona l 



dust, thus neglecting vorticity). In a similar vein, both iRasanenI (120121 ) and iGromed (120111 ) 



assume that time derivatives of metric perturbations ar e small, leading to the conclusion that 



velocity-dependent terms can safely be neglected; and iGromed (120111 ) specifically uses this 
point as a primary argument for dismissing the effects of causal backreaction entirely. Obvi- 
ously, though, if one removes all "causal" aspects from "causal backreaction" , then nothing 
will be left, and it is unsurprising for one to then find that backreaction fails to work, given 
that the principal physics responsible for it has been intentionally set to zero. But this 
basic procedure is conceptually flawed in the case of causal backreaction, which deals with 
propagating modes, since such modes can bring in gravitational perturbation information 
to a local observer from vast spatial volumes, so that a term with a 'small amplitude' can 
actually produce an enormous overall effect. 

The habit of conflating "small" with "negligible" is a continuing error in perturbation 
theory analyses of cosmological evolution, since the size of these effects, in terms of ampli- 
tude, effectively does not matter: it is practically irrelevant how small the time-derivative 
terms may be, when one realizes that such effects are being cumulatively contributed by all 
of the matter within a causal horizon that may be billions of light-years in extent, thus mul- 
tiplying that amplitude by an amount of mass easily large enough to overcome its inherent 
smallness. Furthermore, even if the amplitudes of the relevant perturbative terms for causal 
backreaction were somehow 'magically' made even smaller than they naturally are, this still 
would not shut off the causal backreaction effect, but merely postpone it to begin a little bit 
later, since a somewhat larger causal horizon of self-stabilizing inhomogeneities would then 
be needed to supply a large enough integrated effect to cause an apparent acceleration. In 
that sense, one could never make the amplitudes of those perturbation terms small enough 
to avoid the eventual dominance of causal backreaction, since a big enough causal horizon 
(containing a sufficiently large amount of virializing mass) can always be attained after a 
sufficiently long period of time, which when multiplied by even the smallest amplitude, would 
eventually produce a term of order unity that then proceeds to dominate the cosmic evolu- 
tion. Thus for causal backreaction in an effectively infinite universe, there is no such thing 
as "negligible". 
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3. THE PHENOMENOLOGICAL FORMALISM OF CAUSAL 

BACKREACTION 

We req uire a practical way to implement these ideas, so consider the following. It is 



well known (IWeinberglll972L pp. 474-475) that the Friedmann expansion for some spherical 
volume V can be derived - using nonrelativistic Newtonian equations, in fact - without 
reference to anything outside of that sphere. In a sense, the outside universe might as well 
not even exist for V, in a completely smooth cosmology. 

This changes when perturbations develop, however, since clumped structures external 
to V exert measurable gravitational forces upon the mass within it; and such forces are new, 
as if those clumped structures were brought in from infinity at some recent time during the 
structure formation epoch, to only then begin affecting V. But one cannot feel a gravita- 
tional force, unless one simultaneously feels a gravitational potential; and so the perturbing 
potential within V must also be new in every physically measurable sense as well, approach- 
ing and entering V in causal fashion from some clumped object (and from all others) as they 
develop over time, everywhere in the surrounding universe. 

For a randomly-distributed collection of inhomogeneities, the forces exerted upon the 
mass in V would pull in all directions, largely canceling one another; and indeed, our for- 
malism is based upon a "smoothly-inhomogeneous" universe for which spatial variations are 
essentially averaged out. But what does not cancel out, regardless of how the perturbations 
outside of V are distributed - because it is a non-directional scalar, rather than a vector - is 
their combined contribution to the absolute level of the potential $, which may not matter 
in true Newtonian physics (only potential differences or gradients do), but which does matter 
in general relativity. This overall perturbation potential $sr(^) is treated here as indepen- 
dent of position; but it fully obeys causality, growing as perturbation information comes in 
from the limits of one's observational horizon. Our phenomenological approach thus models 
the inhomogeneity-perturbed evolution of (any volume) V with a metric that contains the 
individual Newtonian-level perturbations from all clumped, virialized structures that have 
been causally 'seen' within V by time t, superposed on top of its (internally-generated) 
background Friedmann expansion. 

This is clearly a very heuristic picture, with the growing perturbation potential $sr(^) 
not representing the literal metric of the universe, but rather acting as a quantifiable short- 
hand for the dominant backreaction effects from structure formation. What is really hap- 
pening, is that virializing overdensities halt their collapse by concentrating their local vor- 
ticity (and/or velocity dispersion); this concentrated vorticity leads to real, extra volume 
expansion, representable (in the final state) at great distances by the tail of a Newtonian 
perturbation potential to the background FRW metric; this Newtonian tail propagates out- 
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ward into space by inducing inward mass flows towards the virialized object from fartlier 
and fartlier distances as time passes; and the total perturbation at any given location will be 
the combined effects of innumerable Newtonian tails of this type, coming in causally (from 
all directions) from within that location's ever-expanding "inhomogeneity horizon". Our 
formalism merely uses $sr(^) to represent all of these complex physical processes in a very 
simple way, and we will see that this is sufficient for producing an observed acceleration fully 
in agreement with observations. 

To phenomenologically model the cosmological phase transition from the extremely 
symmetrical and smooth ( "unclumped" ) state during the early universe, to its almost entirely 
"clumped" state today, we employ a "clumping evolution function", "^{t), defined over the 
range from \E'(t) = (perfectly smooth matter) to \E'(t) = 1 (everything clumped). Note that 
we assume a spatially fiat background - for the pre-perturbed universe - everywhere in our 
calculations, and that we therefore assume the entire universe to be composed of pressureless 
dust with = 1 7^ ^m'^'^ (the superscript "Obs" referring to observational quantities). 

The Newtonian approximation for a single clumped object of mass M (representing a 
virialized, self-stabilized inhomogeneity), embedded at the origin (r = 0) in an expanding, 
spatially fiat , matter-domiii ated (MD) universe can be obtained by linearizing the McVit- 



tie solution (IMcVittid Il933l ). as can be seen from the perturbed FRW expression given in 
Kaloper et al. J2OI0I 1 One gets: 



ds^ ^ -c2[l + (2/c2)$(t)]rft2^[aMD(t)]'[l-(2/c2)$(t)]rfr2 + [aMD(t)]V[rfe2^sin2^rf</>2] , (3) 

where $(t) = {—GM/[aMDit)r]}, and aMD(^) oc t^^^ is the unperturbed MD scale factor. 

Assuming that such mass concentrations will be randomly and (sufficiently) evenly 
distributed throughout the universe, one must integrate over distance (from a particular 
observer's location) for the various clumps, and one must also angle-average over them to 
get the metric for the combined effects that would represent an averaged ds^ typical for that 
observer experiencing displacements in any given direction. Since a displacement of magni- 
tude \dr\ will have a radial component (with respect to any particular mass concentration) 
of \dr\ cos 6*, angle-averaging thus brings in a factor of (cos^ 9) = (1/3) for the spatial part of 
the metric perturbation term, while leaving its temporal part unchanged. 

Consider now the observer to be located at the origin, surrounded by a set of (roughly 
identical) discrete mass concentrations with total mass M, all located within a particular 
spherical shell at coordinate distance r' from the origin, but randomly distributed in {6,(f)). 
From the above arguments, we can write the total, linearly summed and angle-averaged 
metric in 'isotropized' fashion, as: 

ds' = -c2{l - [i?sch(t)/r']} dt' + [aMD(t)]'{l + (l/3)[i?sch(t)/r1} ^rf , (4) 
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where i?sch(t) = {{2GM/c^)/[auD{t)]}, and \dr\^ = {dr^ + r^dO"^ + r"^ sin'^ 9 d^^) = \dx\^ = 
[dx^ + dy"^ + dz^). Regarding the spatial metric perturbation term as a 'true' increase in 
spatial volume, and the perturbation in gu as an 'observational' slowdown of the perceptions 
of observers (relative to the expansion governed by cosmic time t), we see that even if the 
spatial term by itself is not enough to generate a real volumetric acceleration, when coupled 
with the temporal term the entire perturbation may indeed be enough to create a so-called 
"apparent acceleration" sufficient to explain all of the relevant cosmological observations. 

The expression in Equation S] represents the perturbations to the metric due to masses 
at some specific coordinate distance r' - and thus from a specific look-back time t' - as 
seen from some particular observational point. The total metric at that spacetime point 
must be computed via an integration over all possible distances, out to the distance (and 
thus look-back time) at which the universe had been essentially unclustered. Finally, a light 
ray reaching us from its source (e.g., a Type la supernova being used as a standard candle) 
travels to us in a path composed of a collection of such points, where the metric at each point 
must be calculated via its own integration out to its individual "inhomogeneity horizon"; 
and only by calculating the metric at every point in the pathway from the supernova to our 
final location here at r = z = can we figure out the total distance that the light ray has 
been able to travel through the increasingly perturbed metric, given its emission at some 
specific redshift z. 

Consider a light ray emitted by a supernova at cosmic coordinate time t = tsN? which 
then propagates from the supernova at r = tsn, to us at r = 0, t = to- We refer here to the 
geometry depicted in Figure [H 

For each point P = (r, t) of the trajectory, the metric at that point will be perturbed 
away from the background FRW form by all of the virialized clumps that have entered its 
causal horizon by that time. Consider a sphere of (coordinate) radius a, centered around 
point P, with coordinates (a, tret) (where tret < t is the retarded time), defined such that the 
information about the state of the clumping of matter on that sphere at time tret will arrive 
- via causal updating, traveling at the speed of null rays - to point P at the precise time 
t. To compute the fully-perturbed metric at P, we must integrate over the clumping effects 
of all such radii a, from a = out to a^a.^, the farthest distance from P out from which 
clumping information can have causally arrived since the clustering of matter had originally 
begun in cosmic history. 



The full derivations of all formulae shown below are given in iBochnerl (120111 ) ; here we 



quote the resulting equations, which are used for our causal backreaction simulations. 

For a FRW metric with a(t) = aMD(t) = ao{t/tof^^ = c[18tVi^o]^''^ the retarded time 



- 12 - 



tret can be given as: 

t,,t{t,a) = to[{t/toy/^ -af . (5) 

Similarly, we can determine amax? given some initial time tinit at which structure formation 
can be reasonably said to have started (i.e., \E'(t < tmit) = 0): 

amax(t, tinit) = [(t/to)'/' " (tinit Ao) '^'l • (6) 

Now, how the metric is affected at P by a spherical shell of material at coordinate radius 
a depends upon the state of clumping there at the appropriate retarded time: '^[tret(t,a)]. 
The total effect is computed by integrating all shells from a = out to a = «max(^5 ^init); 
but in order to compute the metric perturbation from each shell quantitatively, it is first 
necessary to relate this clumping function to an actual physical density of material. 

Recalling EquationHl we have the perturbation term [Rsch(t)/r'] = {{2GM/c'^)/[r' auDit)]}, 
where the value of M to use here at time t is given by the clumped matter density, [^(t)pcrit(t)], 
times the infinitesimal volume element of the shell. Collecting all relevant factors (and letting 
t' = tj-etit, «)), the integrand at coordinate distance a can eventually be written as: 

[Rsci.it) /r']r>=a-,ia+da) = {12 ^(t') [{to/tf'] [ada]} , (7) 

where for simplification we have used Hq = (2/3)to ^ and the fact that [p{t)a{ty] is constant, 
both true for a matter-dominated universe (but see caveats below). Note also that only \Ef 
is evaluated at the retarded time, tret (^5 since the only "relativistic" piece of propagating 
information which is causally delayed is the state of clumping (\E'[tret(^, «)]) that has just 
then arrived from coordinate distance a to observer P at {r,t). 

From this result, we can now determine the total integrated metric perturbation function 
due to clumping, I(t), as experienced by a null ray (or any observer) passing through point 
P at (r,t): 

/•a!max(t,iinit) 

Iit)= {12 ^[t,etit, a)] [{to/tf'']} ada , (8) 

^0 

with I{t) implicitly being a function of tinn (with I{t) = for t < tmit), as well as of to- 

Finally, we can insert this result back into the formalism of Equation HI to obtain 
the final clumping-perturbed metric that we will use for all of our subsequent cosmological 
calculations: 

ds^ = _c2[i - /(t)] dt' + {[aMD(t)]' [1 + (l/3)/(t)]} \df\' . (9) 

Note that we have implicitly made a significant approximation in this above deriva- 
tion, in that several quantities (e.g., tret, ttmax, [R-Schit) / r']) were calculated relative to the 
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unperturbed background metric, rather than with respect to the inhomogeneity-perturbed 
metric itself. This key approximation is the dropping of what we refer to as "recursive 
nonhnearities" - called that because the integrated function I{t) invariably depends upon 
itself recursively, in an operationally nonlinear way - and this simplification is done here 
out of practical necessity for this initial, proof-of-principle work; but a fully self-consistent 
treatment would need to account for how causal updating is itself slowed by the metric per- 
turbation information carried by causal updating, since the strength of the acceleration as 
2; — )■ should be moderated to some extent by these recursive nonhnearities. Incorporating 
them properly represents a crucial next step for our causal backreaction formalism; but for 
now, we simply note that the output cosmological fits and parameters from our calculations 
here will unavoidably possess some systematic theoretical uncertainty. 

With Equation M as our metric, we must relate unperturbed ("FRW") quantities - 
such as z^^^{t) = [{to/ty^'^ — 1] - to actual physically-observed quantities, such as the 
corresponding observational redshift, given as follows: 



.Obs 



In order to compute the supernova-based luminosity distance function (i.e., the observ- 
able curve which most directly traces out the cosmic expansion history) , we require a formula 
for the coordinate distance r of a supernova going off at coordinate time t, which would be 
seen on Earth precisely now; we find that the resulting expression is: 



.™W^^f{(<«A'P^^^i^}<^<'. (u) 

This formula can then be converted into an expression for the observed luminosity distance 
function, ultimately yielding: 

where Jq = /(to), and t^, 4 dimensionless time ratios (e.g., tr = t/to), with no change 
to the essential form of I{t) (i.e., I{t) = I{tr ■ to) =^ I{tr))- Finally, this luminosity distance 
function di^^p^^t computed for any clumping evolution model \E'(t) is converted to the form 
A(m — M)peTt{z^^^) by subtracting off the c^l function for an empty, coasting universe; and 
the resulting curve can be plotted as a residual Hubble diagram, for use in comparisons 
against theoretical FLRW models and standard candle data from supernova observations. 
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4. MODELS, PARAMETERS, AND RESULTS 

We must design a set of clumping evolution functions, \E'(t^^^), to model the overall 
combined effect due to the linear evolution, the nonlinear regime, and the final virialization of 
individual clustering masses, along with the triggering (often via collisions) of newly-forming 
clumps. 

Observational data and simulations are not sufficiently constraining, so we opt here 
for simplicity. Considering different varieties of time-dependence, we choose: ^(t) oc t^^^, 
proportional to the linear density contrast evolution in a matter-dominated universe (e.g.. 



Kolb fc Turneiill99Cll ): and \E'(t) oc t, the amount of clumping being simply proportional to the 
time that has elapsed. We also try an 'accelerating' clumping rate, \E'(t) o c t^, potentially 



corre sponding to the final nonlinear evolution of a density perturbation flKolb fc Turner 



1990l . p. 322); though the lesser amount of clumping at early times for this model will lead 



to comparatively weak causal backreaction effects. Thus we define: 

to ~ tinit 

^Li„(t) - ^o(^^^) (13b) 

^0 ~ f-init 

^sUt) - ^o(^^)^ (13c) 

where tinit represents the effective beginning of clumping (i.e., \E'(t < tinit) = 0), and \E'o 
represents the ultimate degree of clumping today. 

Appropriate values for the two adjustable parameters, z-mit (= [(^o/^imt)^^'^ — 1]) and "^o, 
are chosen via astrophysical considerati ons. First, we link the beginning of structure forma- 



tion to cosmological reionization. Since iDunkley et al.l ( 120091 ) supports a possibly extended 



period of partial reionization lasting from 2; ~ 20 — 6, we conservatively bracket this by using 
^^init = (5, 10, 15, 25) for our simulation runs. 

To specify \E'o, we take typical concordance values such as f^M^*^ = 1 — Q'^^^ ~ 0.27, with 
Qf^^ ~ 0.04 (thus f^DM ~ 0.23), and scale them up to flatness without Dark Energy. This 
yields Qf^"^ ~ [0.04(1.0/0.27)] ~ 0.15, and ~ 1-^]^^^ ^ o.85. Estimating the various 

fractions of da. rk and baryonic matter which may currently be clumped versus unclumped 



( lBochnerll201l[ l. we choose the range \E'o = (0.78, 0.85, 0.92, 0.96, 1.0) for our runs. 



Four different Zi^it values, five different \E'o values, and three different clumping evolution 
models gives us 4 x 5 x 3 = 60 simulation runs in total. Residual Hubble diagrams of 
these model cosmologies show them to be broadly successful in reproducing a Cosmological- 
Constant-like observed acceleration, with the \E'sqi. runs resembling flat ACDM with Q\ ~ 
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0.3 — 0.4, the \1/Lin runs resembling ~ 0.5 — 0.8, and the \1/md runs resembhng ~ 
0.65 — 0.97. Quahtatively, a ~dozen of these model cosmologies fit the data essentially as 
well as the best-fit fiat ACDM model. Two panels containing several of these best runs are 
shown in Figures |2] and [31 

For our quantitative analysis below, we fit the residual distance modulus function for 



each model, Afipcn = A(m — M)pert, to the SCP Union compilation (IKowalski et al.l 120081 ) 



of 307 Type la supernovae (SNe). Using their publicly available /isN,i and csN.i data, we 
compute values of Xfh these A/zpert curves, minimized for each run by optimizing H^^^. 
Given these results, the corresponding fit probability values (-Ppit) for each run are then 
computed using the distribution with A^doF = (307 — N^n) degrees of freedom, where 
A^Fit = 3 for our ^(t) models (optimizing ^o, Zi^it, and if,?'^"), A^pit = 2 for fiat ACDM 
(optimizing Q;^ = 1-Qm and i^o^^'), and A'pit = 1 for fiat SCDM (optimizing H^^"" only). 



FORGING A NEW CONCORDANCE 



Beyond just explaining the apparent acceleration seen in the SNe data, a true model of 
the universe must satisfy a variety of complementary constraints before it can be regarded 
as a fully-consistent replacement for the traditional ACDM "Cosmic Concordance" . 

To obtain a new concordance without Dark Energy, we must relate the observational 
("dressed") parameters to the unperturbed FRW ( "bare" ) parameters of our models. These 
mathematical relationships are derived in iBochnerl ( 1201 ll ): here we will simply quote the 
output results for our best-fitting models. 



rFRW 



versus measured Hubble Constant Hq^^); the 



The parameters studied include: H^ 
physical age of the universe, t^^*^; the unperturbed cosmic density, nff^; the observed CMB 
acoustic scale, l^^^; the apparent (overall) cosmic equation of state, Wq^^; and the observed 
jerk (or jolt) parameter, j^^^. Observables H^^^, Wq^^, and Jq^^ are computed from the 
Taylor series expai ision coefficients of the luminosity distance function, di^, pertiz'^^^), for 



,Obs 



(see, e.g.. lRiess et al.l 120041 ). 



Besides achieving an 'acceleration 
batio ns exist (i.e., I(t) > ), one has H^ 
(e.g. 



FRW 




< 



1/3), our key result is that when pertur- 
which is astrophysically interesting 



rrObs 



Blanchard et al.ll2003[ ) and solves a number of problems with matter-only cosmologies. 



First, since (t™^, t?^^) oc 1/if™^, a lo w if™^ can solve the classic Age Problem/Crisis 



in cosmology (e.g.. lKolb fc Turnerlll990l ) by increasing t, 



^Obs fj^Qj^ ^]^g 



from always-decelerating flat SCDM cosmologies, to values of ' 
perturbed models. 



9-10 GYr expected 
13 — 14 GYr for our strongly- 
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Secondly, since = pu/ Pcrit oc pm/Hq, a universe that appears to have insufficient 
density for closure due to oc pm/ {H^^^Y ~ 0.3, may very well be spatially flat (in a pre- 
perturbed sense) for the same physical matter density pm, with oc PmJJHo^Y — 1- 



(And note that CMB- related tests of cosmic flatness - e.g., iLarson et al.l (1201 ll ) - primarily 
measure the pre-perturbed era.) Thus we can reconcile the apparent contradiction between 
the CMB measurements of l^^^ indicating flatness, and our seemingly low-matter-density 
universe (as indicated by the growth of structure), without requiring the addition of any 
non-clustering Dark Energy species to ffil the apparent gap between ^l^^ ~ 0.3 and fixot = 1- 

Examining our 60 simulated cosmological models quantitatively, we informally choose 
a set of 'best' runs - six "^un models and six \I/md models - which provide very good SNe 
data fits, while simultaneously producing good cosmological parameters. A truly optimized 
search over the (\l/o, ^^mit) parameter space is not really called for at this toy-model stage of 
our formalism; but a quick effort at optimization reveals an extensive 'trench' in parameter 
space for the \1/md runs with quite low (and highly degenerate) Xpit values, yet offering 
a wide variation to choose from regarding their output cosmological parameters. So for 
illustration purposes, we arbitrarily select one low-x|;^ case, (\l/o; -^^init) = (0.768,14), as a 
so-called "semi-optimized" run, giving us a total of thirteen "best runs" for more in-depth 
study. 

Residual Hubble diagrams of these best runs are shown in Figure HI These thirteen 
cosmological models clearly produce good Hubble curves, being visually almost indistin- 
guishable from one another (and from the now best-fit, Q/^ = 0.713, Concordance ACDM 
model) in the SN-data-rich region of 2°^^ ~ 0.1 — 1. The causal backreaction formalism 
has therefore succeeded in reproducing the apparent cosmic acceleration as it is actually 
measured, via SNe standard candles. 

More quantitatively, the comprehensive output data from these thirteen best runs (and 
from best-fit ACDM and SCDM, for comparison) are given in Table [H 

Summarizing the results: these runs fit the SNe data essentially as well as ACDM in 
terms of Xpif comparable in Ppit- The mid-range values of Jq indicate a strong 

enough backreaction effect to explain the apparent acceleration, without being so large (i.e., 
Jo ~ 1) as to indicate a substantial breakdown of our approximation of gravitational linearity 
in simply summing together individual metric perturbations to produce the total overall 
perturbation. 

Crucially, these models all solve the Age Problem without a Dark-Energy-induced ac- 
celeration, with t^^^ ~ 13.2 — 14.5 GYr; and they also achieve (pre-perturbation) spatial 
flatness to within a reasonable degree of uncertainty. Specifically, the "^un runs yield 
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^kIP^ ~ 0.90 — 1.04, very close to unity; and while the fi^^^ values for several of the 
\E'md runs are somewhat high, lowering the (pre-adopted) normalization of = 0.27 to 
the (also reasonable) value of = 0.24 (and dropping the worst ^E'md run) moves them 
down to f^M^^ ~ 0-88 ~ 1.25. These ~dozen runs also succeed fairly well in bracketing the 
SNe-best-fit ACDM value (for no radiation) of ^aIacdm = 285.4 for the CMB, performing 
vastly better than anything achievable with a matter-only open (i.e., ~ 0.3) universe. 

Possessing very negative effective equation of state values, Wq^^ < —0.7, these models 
clearly produce a strong amount of apparent acceleration, despite the absence of a negative- 
pressure cosmic component; and we note that since t he precise cosmological parameters 



obtained from SNe fits are highly model-dependent (e.g.. ICattoen &: Visserll2008l ). it is there- 
fore unnecessary for alternative-cosmology models to exactly reproduce the w^^^ result from 
Concordance ACDM (e.g., Wq^^ = —0.713), but merely to generate an 'accelerative' value of 
Wq^^ strong enough to provide a sufficiently good fit to the supernova data (as has indeed 
been done here). Taken together, these results of a sufficiently-strong apparent acceleration, 
quantitatively good SNe fits, and good cosmological parameters all indicate that we have 
found a selection of causal backreaction models that successfully achieves the basic goals of 
an alternative cosmic concordance, without the use of any form of Dark Energy. 

Lastly, consider the observed jerk parameter, Jq^^. Since fiat ACDM cosmologies pos- 
sessing only dust and a Cosmological Constant have j{z) = jo = 1 for all time (for any choice 
of fl\ = 1 — ^m), searching for deviations from jo = 1 represents a fundamental test of the 
FLRW-ACDM paradigm, itself. Coming from the third-order series term for di^{z^^^), j^^^ 
is not yet observationally well-constrained, and so can be considered a testable prediction of 
our formalism. 

From Tabled] we have j^^^ ~ 2.5 — 5.5 for our thirteen best runs, a range that can be 
further narrowed down to j^^^ ~ 2.6 — 3.8 by restricting consideration to the six of those 
runs which are most consistent with all other cosmological parameter constraints (specifically: 
^Lin with (^0, ^init) = (1-0, 25), (1.0, 15), and (0.96, 25); and ^md with (^o, ^init) = (0.78, 10), 
(0.85,5), and (0.92,5)). These jerk parameter values are clearly in excess of the ACDM 
prediction of jo = making a comparison with observational results imperative. Though 
the great difficulty of precisely measuring this parameter has led to a wildly- v a rying range 



of estimates, one potentially reliable estimate has been adapted in iBochnerl ( 1201 ll ) from 



the detailed results (D. Rubin, 2010, private communicati on) of the (SNe-|-CMB- |-BAO) 



analysis done with the SCP Union2 SNe compilation data (lAmanullah et al.l 120101 ). using 
statistical and systematic errors, which in terms of the jerk parameter yields the result: 
j^^" ± Aj^^" ~ 2.129 ± 1.292; i.e., j^^' ~ 0.8 - 3.4. While still consistent with Concordance 
ACDM, it is focused predominantly on the high side, as generally predicted by our formalism 
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(if not quite so high as several of our models predict); an observationally interesting situation, 
being not yet conclusive in either direction. 

One important caveat, however, is that we cannot conclude with certainty - based 
only upon the results in this paper - that high jerk parameters are a robust feature of our 
formalism. Too many theoretical uncertainties exist within this toy-model approach to be 
sure, just yet, that this is an 'iron-clad' prediction of causal backreaction. This includes our 
neglect here of both recursive and gravitational nonlinearities; as well as the oversimplified 
nature of these "^(t) clumping evolution models, which neglect astroph ysical processes such a s 



feedback from star formation and the shock- heating of baryons (e.g., ICen &: Ostrikerll2006l ). 
All such effects (except perhaps gravitational nonlinearities) should result in the damping of 
the apparent acceleration, especially at later times, resulting in Jq^^ values that are reduced 
compared to those seen in the output results of the naive models presented here. Thus 
the theoretical situation is as uncertain as the observational situation is with regard to this 
key cosmological parameter, and more theoretical development of the causal backreaction 
formalism needs to be done before this paradigm can be conclusively tested by measurements 
of Jq^^. It is possible, in any case, that other means may be found to test it; for example, 
the fact that causal backreaction does not occur in completely smoothly-distributed fashion, 
but is instead most concentrated near virializing masses, may lead to interesting feedback 
behaviors affecting the formation of stars, galaxies and galaxy clus ters. Conceivably, this 



may have some useful application to issues such as galaxy downsizing fICowie et al.lll996r). the 



cuspy CDM halo problem, and the possible dearth of dwarf satellite galaxies (e.g., iPrimack 



20031 ): but such connections are speculative, and finding a variety of concrete methods for 



testing causal backreaction remains of great importance for this paradigm. 



6. SUMMARY AND CONCLUSIONS 

Our starting point in this paper has been the search for an alternative to Dark Energy for 
attaining a Cosmic Concordance, without resorting to the added complexity of modifications 
to gravity, or the non-Copernican 'specialness' implied by a local cosmic void. The problems 
with Dark Energy are well known, including the magnitude and coincidence fine-tuning 
problems for a pure-A Cosmological Constant, and the stability problem for a more dynamical 
form of Dark Energy possessing the ability to cluster spatially. The latter (DDE) case 
requires the ad-hoc addition of some form of new pressure term to support it against collapse, 
given the locally attractive nature of the negative pressure required by the DDE to power the 
cosmic acceleration. Such a new pressure term, if adiabatic (e.g., degeneracy pressure), would 
represent a form of positive pressure contributing to a cosmic deceleration that partially or 



- 19 - 



totally nullifies the acceleration meant to come from the DDE; and if non-adiabatic (due to 
some new effects from the DDE Lagrangian) would invalidate the Dark Energy as a perfect 
fluid, thus calling into question the usual cosmic implications of w^^ < in the FLRW 
acceleration equation. 

Emphasizing the importance of the fact that 'negative pressure' is locally attractive in 
character - rather than repulsive, as it is often regarded and popularly described - and that 
normal gravitational attraction therefore represents (at least in a non-technical sense) a form 
of negative pressure, we have therefore made a virtue of necessity by recruiting cosmological 
structure formation itself, based upon ordinary gravitational forces, as the driver of the 
observed (possibly apparent) acceleration. This approach, known generally in the literature 
as "backreaction" , removes all of the aforementioned problems by eliminating the need for 
any Dark Energy species, while solving the Coincidence Problem by naturally linking the 
onset of cosmic acceleration to the emergence of cosmic structure. This linkage inevitably 
leads to the creation of observers (such as ourselves) just in time to see this 'coincidence'. 

Noting that different varieties of backreaction proposed by other authors have been 
largely unsuccessful in their attempts to account for the observed acceleration, and that 
powerful arguments have been advanced which claim that backreaction as a paradigm can- 
not be made strong enough to succeed at this task, we have explained why each of these 
"no-go" arguments arc functionally invalid, due to a number of astrophysically inappropriate 
approximations. These include the dropping of (squared) vorticity and velocity dispersion 
terms in backreaction calculations; the unrealistic boundary conditions and isolation of over- 
densities in Swiss-Cheese models; and the neglect of Newtonian-level perturbations, due to 
well-known results from an influential backreaction formalism (and from perturbation the- 
ory in general) which adopt an inappropriate 'slow-motion' approximation, and thus lack 
essential effects due to the causal propagation of gravitational information from developing 
inhomogeneities ( "causal updating" ) . 

Informed by the consequential but flawed assumptions of such no-go arguments, we have 
developed a phenomenological formalism of "causal backreaction" , in which the apparent ac- 
celeration of the universe is the result of the universal phase transition from smoothness to 
dumpiness, and is achieved through the combined perturbative contributions of all viri- 
alizing, self-stabilizing cosmic structures to the metric of a typical cosmological observer. 
Averaging over location and direction, the evolution of this "smoothly-inhomogeneous" cos- 
mology is determined by a single input function, \E'(t), chosen from a selection of clumping 
evolution models which are specified heuristically but guided by observed astrophysical data. 

Using a set of such clumping evolution models with physically-motivated input param- 
eters, we have conducted 60 simulation runs, finding a ~dozen-plus solutions that fit the 
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SCP Union Compilation Type la supernova data essentially as well as the best-fit ACDM 
model. Furthermore, about half of these models give good values (within reasonable the- 
oretical and observational uncertainties) for all other calculated cosmological parameters 
studied here, such as: the age of the universe (thus solving the Cosmic Age Problem); the 
matter density as a fraction of the unperturbed FRW critical density (thereby achieving 
effective spatial flatness with matter alone); the deceleration parameter (therefore achieving 
a sufficiently-strong apparent acceleration) ; and the characteristic angular scale of the CMB 
acoustic peaks (thus matching observations from the early, pre-structure-formation cosmic 
epoch). All of these goals are achieved without altering the measured physical density of 
cosmic matter, and without the introduction of any form of Dark Energy, thus forming the 
basis of an alternative concordance for a matter-only universe. 

Looking for testable ways to distinguish our causal backreaction formalism from Cosmo- 
logical Constant ACDM (and from other slowly-evolving forms of Dark Energy not too far 
from A), we have seen that while the jerk parameter must always obey the condition j^^^ = 1 
for flat ACDM, the results of those of our causal backreaction models which successfully fit 
the supernova data instead predict j^^^^ values strongly in excess of unity, a result which 
has some support from the observational data. Nevertheless, we note these caveats: ACDM 
still remains consistent at the la level with most observations; the definitive observational 
constraints on j^^^ remain extremely weak, in any case; and perhaps most importantly, 
theoretical uncertainties for the highly simplified causal backreaction formalism introduced 
here has likely led to overestimates in our predicted values of the jerk parameter, a potential 
problem which must be remedied through the further theoretical development of the model. 

Expressing the conclusions of this paper in one sentence: if someone were to ask, "What 
is the force behind the cosmic acceleration?", the answer we would give is that it is not a 
'force' at all; rather, the motivating factor is the total, summed effect of the Newtonian 
tails of individual metric perturbations, produced by the virialization of innumerable self- 
stabilizing structures filling the universe, with these influences propagating causally towards 
all observers from the extreme edges of their observable cosmic horizons. 

I am grateful to Jacob Bekenstein, Krzysztof Bolejko, Varoujan Gorjian, Wayne Hu, 
Marek Kowalski, Edvard Mortsell, David Rubin, Ran Sivron, David Wiltshire, and Ned 
Wright for brief but helpful communications; and I am especially grateful to Arthur Lue for 
several helpful and clarifying discussions. 
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(r,t) = (0,Q 



Fig. 1. — Geometry for computing the inhomogeneity-perturbed metric at each point along 
the integrated path of a hght ray from a supernova to our observation point at Earth. 




Fig. 2. — Residual Hubble diagrams for selected \1/Lm(^) clumping functions (solid lines), 
plotted versus fiat ACDM {^a = 0.73 = 1 - ^m) and SCDM (^m = 1) cosmologies. 
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Fig. 3. — Residual Hubble diagrams for selected ^md(^) clumping functions (solid lines), 
plotted versus flat ACDM (^a = 0.73 = 1 - ^m) and SCDM (Qm = 1) cosmologies. 
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Fig. 4. — Residual Hubble diagrams for our thirteen 'best' models, shown along with the 
best-fit flat SCDM and Concordance ACDM (Qa = 0.713) cosmologies, and the binned and 
averaged SCP Union SNe data. From highest to lowest (at z^^^ ~ 10), the plotted *md(^) 
curves have the parameters: (*o,-2init) = (0.96,5), (0.92,5), (0.85,5), (0.78,10), (0.78,15), 
(0.768, 14), and (0.78,25). From highest to lowest (at z^^"" ~ 10), the plotted ^Lin(t) curves 
have the parameters: (^o.^init) = (1-0,10), (1.0,15), (0.96,15), (1.0.25), (0.96,25), and 
(0.92,25). The individually-optimized H^^^ values for each model are depicted as differing 
vertical offsets for the curves. 
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Tal)lo 1: Cosiiiolog 


ic-al Paraiiieters Deri\XHl from 


our 'Best' Siinulatiori 


Runs 






(*0, -^init) 


Xlit 


p a 


lo 


^.Obs b 


^Obs c 


^FRW d 


j-Obs e 


OFRW f 


Obs 


,Obs 

Jo 


/Obs 
'A 


^Lin Clumping Model Runs 


(1.0,25) 


312.1 


0.362 


0.57 


1.12 


69.95 


40.68 


13.56 


1.037 


-0.817 


3.45 


284.2 


(1.0,15) 


313.3 


0.344 


0.55 


1.12 


69.64 


41.63 


13.45 


0.971 


-0.784 


3.11 


287.6 


(1.0,10) 


315.4 


0.314 


0.52 


1.12 


69.27 


42.80 


13.30 


0.897 


-0.746 


2.76 


291.3 


(0.96,25) 


314.8 


0.323 


0.55 


1.11 


69.38 


41.57 


13.38 


0.968 


-0.759 


2.92 


285.9 


(0.96,15) 


316.6 


0.297 


0.52 


1.11 


69.11 


42.46 


13.28 


0.911 


-0.732 


2.67 


289.1 


(0.92,25) 


319.0 


0.265 


0.53 


1.11 


68.85 


42.44 


13.21 


0.905 


-0.706 


2.49 


287.6 


^'md Clumping Model Runs 


(0.78,10) 


312.1 


0.362 


0.63 


1.12 


69.96 


38.25 


13.88 


1.204 


-0.801 


3.15 


277.5 


(0.78,15) 


312.2 


0.360 


0.68 


1.12 


70.84 


36.14 


14.17 


1.409 


-0.895 


4.17 


270.5 


(0.78,25) 


316.8 


0.295 


0.72 


1.12 


71.80 


34.22 


14.46 


1.642 


-1.001 


5.51 


263.7 


(0.85,5) 


313.9 


0.336 


0.56 


1.13 


69.48 


41.21 


13.60 


0.991 


-0.747 


2.59 


288.7 


(0.92,5) 


312.1 


0.363 


0.60 


1.14 


70.71 


39.41 


14.00 


1.144 


-0.871 


3.75 


285.6 


(0.96,5) 


315.5 


0.313 


0.63 


1.14 


71.52 


38.36 


14.25 


1.248 


-0.954 


4.67 


283.8 


Semi- Optimized ^md Clumping Model Run 


(0.768,14) 


311.7 


0.369 


0.66 


1.12 


70.37 


36.91 


14.03 


1.324 


-0.845 


3.63 


272.4 




Comparison Values from Best-Fit flat ACDM Model {Qf^ 


= 0.713 


= 1-11 


m) 






311.9 


0.380 




1.0 


69.96 


69.96 


13.64 


0.287 


-0.713 


1.0 


285.4 


Comparison Values from Best-Fit flat SCDM Model (Qa = 0, = 1) 




608.2 


3.4E-22 




1.0 


61.35 


61.35 


10.62 


1.0 


0.0 


1.0 


287.3 



"Ppit values are calculated using A^doF = 304 for the ^'Lin and ^'md models, A^doF = 305 for flat ACDM, and 
A^DoF = 306 for flat SCDM. 

''Each z'^^^ value quoted here corresponds to z^^^ = 1. 
°A11 Hq values are given in km s~^Mpc~^. 

''Each 7?™^ is computed relative to the corresponding optimized Hq^^ value for that run. 
"^tg values are listed in GYr, and all quantities are computed assuming no radiation (Qr = 0). 
•'^AU O^^^ values given for the '^un and ^'md models are normalized to Sl^^ = 0.27. 



